arXiv:hep-lat/0112027vl 17 Dec 2001 


1 


Quark propagator from an improved staggered action in Laplacian and 
Landau gauges 

Patrick O. Bowman’^*, Urs M. Heller*^ and Anthony G. Williams'^ 

'^Department of Physics and CSIT, Florida State University, Tallahassee FL 32306-4120, USA 
*’CSSM and Department of Physics and Mathematical Physics, Adelaide University, Australia 5005 

Studies of gauge dependent quantities are afflicted with Gribov copies, but Laplacian gauge fixing provides 
one possibie solution to this problem. We present results for the lattice quark propagator in both Landau and 
Laplacian gauges using standard and improved staggered quark actions. The standard Kogut-Susskind action 
has errors of 0{a^) while the improved “Asqtad” action has C9(a^), 0{a^g^) errors and this improvement is seen 
in the quark propagator. We demonstrate the application of tree-level corrections to these actions and see that 
Landau and Laplacian gauges produce very similar results. In addition, we test an ansatz for the quark mass 
function, with promising results. In the chiral limit, the infrared quark mass, M{q^ = 0) is found to be 260 db 20 
MeV. 


1. The quark propagator 

The quark propagator is a fundamental quan¬ 
tity of QCD. Though gauge dependent, it mani¬ 
festly displays dynamical chiral symmetry break¬ 
ing, contains the chiral condensate and Aqcd, 
and has been used to compute the running quark 
mass Some model hadron calculations rely 
on ansatze for the quark propagator, yet on the 
lattice we have the opportunity to study it in a di¬ 
rect, nonperturbative fashion. Quark propagator 
studies can be complicated, however, by strong 
lattice artefacts ||,^; it has also been studied us¬ 
ing the overlap action [|j . 

We are required to fix a gauge and we choose 
the Landau and the Laplacian gauges . We use 
Wilson glue at /3 = 5.85 (a ~ 0.125 fm) on a 16^ x 
32 lattice and six quark masses from am = 0.075 
down to 0.0125 (115 to 19 MeV). Calculations 
were done on 80 configurations. 

We use the standard Kogut-Susskind (KS) ac¬ 
tion and the “Asqtad” quark action , a fat-link 
staggered action that combines three-link, five- 
link and seven-link staples to minimise flavour 
changing interactions along with the three-link 
Naik term and planar five-link Lepage term. The 
coefficients are tadpole improved and chosen to 


‘Presented by POB. 


remove all tree-level 0{a^) errors. This action 
was motivated by the desire to improve flavour 
symmetry, but has also been reported to have 
good rotational properties and small mass renor¬ 
malisation 0- 

In the (Euclidean) continuum, Lorentz invari¬ 
ance allows us to decompose the full quark prop¬ 
agator into Dirac vector and scalar pieces 

S-^{p^) = Z-^{p^)[i^-p + M{p^)]. ( 1 ) 

Asymptotic freedom means that, as p^ —> oo, 
—>■ ij-p + m, (the free propagator) where 
m is the bare quark mass. 

From consideration of the tree-level forms of 
our two lattice actions, we define the momentum 
variables = sin(p^) for the KS action and 

= sin(p^) [1 + ^ sin2(T/x)] (2) 

for the Asqtad action, where is the usual lattice 
momentum. 
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By considering the propagator as a function of 
instead of p^, we ensure that the lattice quark 
propagator has the correct tree-level form and 
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hopefully better approximates its continuum be¬ 
haviour. This is the same philosphy that has 
been used in studies of the gluon propagator (see 
Ref. and references therein). See also foot¬ 
note 6 in Ref. Q. 

To help us identify lattice artefacts - as we are 
employing only one set of configurations - we sep¬ 
arate the data according to the direction in which 
the momentum lies. Data from momenta lying 
wholly on a spatial cartesian direction are plot¬ 
ted with squares, along the temporal direction, 
triangles, and along the four-diagonal, with dia¬ 
monds. Circles represent any other combination. 

2. Laplacian gauge 

Laplacian gauge is a nonlinear gauge fixing that 
respects rotational invariance yet is free of Gribov 
ambiguities. Although it is difficult to understand 
perturbatively, it is equivalent to Landau gauge 
in the asymptotic region It is also compu¬ 

tationally cheaper then Landau gauge. There is, 
however, more than one way of obtaining such a 
gauge fixing in SU{N). The three implementa¬ 
tions of Laplacian gauge fixing used here are: 

1. 9^(1) gauge (QR decomposition), used in 
Ref. 

2. 9^(11) gauge, where the complex 3x3 matrix 
is projected onto SU{i) by maximising its 
trace. This will be discussed in more detail 
in an upcoming paper on the gluon propa¬ 
gator Q. 

3. 9^(111) gauge (Polar decomposition), the 
original prescription described in Ref. . 

The gauge transformations employed in Lapla¬ 
cian gauge fixing are constructed from the low¬ 
est eigenvectors of the covariant lattice Laplacian 
operator. The way that these eigenvectors trans¬ 
form under gauge transformations provides the 
uniqueness of the Laplacian gauge. The three im¬ 
plementations discussed differ in the way that the 
gauge transformation is constructed from those 
eigenvectors. We can think of each projection 
method as defining its own Laplacian gauge. In 
all cases the resulting gauge is unambiguous for 
all configurations except a set of measure zero. 
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Figure 1. The quark Z function for the Asqtad 
action (ma = 0.05). When we use the momentum 
q, defined by the tree-level form of the action, 
rotational symmetry breaking is reduced. 

3. Results 

3.1. Performance of the Asqtad action 

We investigated the applicaton of tree-level cor¬ 
rection to the quark propagator by comparing the 
Z functions using p and q. One example is shown 
in Fig. 1^. With the Asqtad action we see that 
hypercubic artefacts are small in any case, but 
when we use q instead of p they nearly vanish. 

It is less clear which momentum variable should 
be used for the mass function, because at tree- 
level it is not multiplied by the momentum, but 
for consistency we use q here as well. In the case 
of the mass function, the choice of momentum 
will actually make little difference to our results. 

In Fig. H the mass function is plotted, in Lan¬ 
dau gauge, for both actions with quark mass 
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Figure 2. Mass function for quark mass ma = 
0.05 (m ~ 77 MeV), KS action (top) and Asqtad 
action (bottom) in Landau gauge. 


ma = 0.05. We see that the KS action gives 
a much larger value for M(0) than the Asqtad 
action and is slower to approach asymptotic be¬ 
haviour. The Asqtad action also shows slightly 
better rotational symmetry. 

The Asqtad action displays clearly better ro¬ 
tational symmetry in the quark Z function and, 
curiously, improved infrared behaviour as well. 
The Asqtad action also displays a better approach 
to asymptopia, more smoothly approaching one 
in the ultraviolet. Furthermore, the relative im¬ 
provement increases as the quark mass decreases. 
Comparing the mass function for the two actions 
at ma = 0.0125, the lowest mass studied here, the 
propagator has much less infrared noise with the 
Asqtad action than with the KS action. 
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Figure 3. Gauge dependence of the quark Z (top) 
and mass (bottom) functions for the Asqtad ac¬ 
tion (ma = 0.05). Points marked with open cir¬ 
cles are in Landau gauge and solid triangles are 
in 9^(1) gauge. Data has been cylinder cut. 

3.2. Comparison of the gauges 

Fig. ^ (top) shows the Z function for the Asq¬ 
tad action in Landau and 9^(1) gauges. Data has 
been cylinder cut | pd| for easier comparison. They 
are in excellent agreement in the ultraviolet but 
differ significantly in the infrared. The dipping 
of Z in the infrared is associated with dynamical 
chiral symmetry breaking. There appears to be 
some slight difference in the Z function between 
9^(1) and 9^(11) gauges. 

Fig. I (bottom) shows the mass function for 
the Asqtad action in Landau and 9^(1) gauges. 
The two mass functions agree in the ultraviolet 
and in the infrared, but the Landau gauge mass 
function sits slightly higher in the intermediate 
region. The mass functions are nearly identical in 
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Figure 4. Comparison of the quark Z functions 
for the three quark masses ma = 0.0125, 0.025 
and 0.05, with the Asqtad action in 9^(11) gauge. 
Data has been cylinder cut. As in Landau gauge, 
they agree to within errors, although there is a 
systematic ordering of the infrared points from 
heaviest quark (top) to lightest (bottom). 

9^ (I) and 9^ (II) gauges. We have also found that 
in Landau gauge the mass function has slightly 
less anisotropy at this lattice spacing. 

Landau gauge seems to respond somewhat bet¬ 
ter than 9^(11) gauge to vanishing quark mass. 
With the smallest quark mass (ma = 0.0125) 
the lowest momentum points have large errors in 
9^(11) gauge compared with Landau rauge. This 
can be seen in the Z function in Fig. jj. 

9^(111) performs very poorly. We found that 
many of the matrices had vanishingly small de¬ 
terminants (compared to numerical precision), 
which destroyed the projection onto SU{3). 
Problems with 9^(111) have also been seen in the 
gluon propagator fl^ . 

3.3. Modelling the quark propagator 

We performed an extrapolation to the chiral 
limit using a quadratic fit. The model ansatz 

+ (4) 

which is a generalisation of the one used in 
Ref. was then fit to each mass function. The 
resulting fit parameters are shown in Table and 
Fig. I shows the best fit for the mass function in 
the chiral limit. We show here only results for 
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Figure 5. Mass function extrapolated to the chiral 
limit. Errors are Jack-knife. Fit parameters are 
c = 0.030(4), A = 870(60) MeV, mo = 0.0, a = 
1.52(23), / dof = 0.49. 

Landau gauge, as they are compatible with the 
Laplacian gauge results, but slightly cleaner. We 
find that a > 1 is increasingly favoured as the 
quark mass approaches zero. 

4. Conclusions 

We have seen that the lattice quark propaga¬ 
tor has better rotational symmetry and displays 
more rapid approach to asymptotic behaviour 
with the Asqtad action than with the standard 
Kogut-Susskind action. Three implementations 
of the Laplacian gauge were investigated, and 
it was found that 9^(1) and 9^(11) gauges gave 
similar results to Landau gauge. 9^(111) worked 
very poorly. The mass function showed very lit¬ 
tle sensitivity to the choice of gauge, but some 
change was seen in the quark Z function. We at¬ 
tempted to model the mass function and saw that 
the ansatz provided a good fit to the data. 

As we have simulated on only one lattice, it re¬ 
mains to do a thorough examination of discreti¬ 
sation and finite volume effects. A natural exten¬ 
sion of the mass function ansatz would be to in¬ 
clude the correct asymptotic behaviour, and this 
would require testing on much finer lattices. 
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Table 1 

Best-fit parameters for the ansatz, Eq. (^), in Landau gauge, in physical units. The table is divided into 
two sections, the first has a = 1.0 fixed and the second leaves a as a free parameter. For the last fit in 
each case, the ultraviolet mass, mg, was fixed to zero. Generally, a is not well determined by the fits, 
but a > 1 seems to be favoured in the chiral limit. 


m (MeV) 

c 

A (MeV) 

Too (MeV) 

a 

M(0) (MeV) 

/dof 

115 

0.40(2) 

920(20) 

144(7) 

1.0 

467(9) 

0.38 

96 

0.36(5) 

890(70) 

118(8) 

1.0 

440(20) 

0.42 

77 

0.41(5) 

830(70) 

95(7) 

1.0 

430(20) 

0.42 

58 

0.45(4) 

770(50) 

70(7) 

1.0 

420(20) 

0.51 

38 

0.49(6) 

720(60) 

44(6) 

1.0 

400(30) 

0.56 

19 

0.54(7) 

670(60) 

18(6) 

1.0 

380(30) 

0.69 

0 

0.56(8) 

660(50) 

-12(6) 

1.0 

350(30) 

0.66 

0 

0.77(15) 

530(50) 

0.0 

1.0 

410(40) 

1.3 

115 

0.28(1) 

1000(30) 

157(7) 

1.25(4) 

433(7) 

0.38 

96 

0.28(2) 

975(40) 

131(9) 

1.26(10) 

408(9) 

0.37 

77 

0.30(5) 

940(50) 

110(10) 

1.29(6) 

380(10) 

0.36 

58 

0.30(2) 

920(40) 

80(6) 

1.30(2) 

360(10) 

0.41 

38 

0.36(6) 

830(50) 

55(7) 

1.28(7) 

350(20) 

0.46 

19 

0.34(4) 

830(100) 

28(5) 

1.35(14) 

310(40) 

0.55 

0 

0.30(4) 

870(60) 

0.0 

1.52(23) 

260(20) 

0.49 
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